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Quantum Heisenberg spin chains with random couplings and spin sizes are studied using a 
real-space renormalization group technique. These systems belong to a new universality class of 
disordered quantum spin systems realized in e.g. SrsCuPti-^Ir^Oe- The low-energy fixed point is 
characterized by the formation of weakly coupled large effective spins. At low temperature T the 
entropy obeys a power law T 6 (S ta 0.44), and the susceptibility follows a Curie-like behavior. 
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A class of quasi-one-dimensional (ID) spin systems 
with an interesting type of disorder was recently discov- 
ered in the compound Sr3CuPti_ x Ir x 06. In the pure 
material, the Cu (spin S — |) and Pt (S — 0) al- 
ternate along chains and the Cu spins interact antifer- 
romagnetically with each other |IJ. In the compound 
Sr3CuPti_ x Ir x 06 a fraction of the spinless Pt atoms is 
replaced by Ir with a spin S = | that couples ferro- 
magnetically to the neighboring Cu spins. Thus this al- 
loy represents a quasi-lD system with random ferromag- 
netic (FM) and antiferromagnetic (AF) bonds of fixed 
strength |2|]. Measurements of the uniform susceptibil- 
ity as a function of temperature (2K < T < 300K) 
revealed a crossover around 10K between the standard 
high-temperature and a low-temperature Curie behavior 

In a recent theoretical study we considered a related 
model of a quantum spin chain with random FM and 
AF nearest-neighbor coupling of equal strength J [|J. 
The analysis of high-temperature series for the suscepti- 
bility compares qualitatively well with the experimental 
data H. Our analysis of the entropy led to the sur- 
prising result that a considerable fraction of the degrees 
of freedom remains uncorrelated down to a temperature 
ksT ~ J/5. This unusual property together with the 
Curie-like behavior of the susceptibility indicated the ex- 
istence of very weakly coupled effective spin degrees of 
freedom in the low-energy regime, which we identify in 
the following way. Our system consists of a sequence of 
alternating FM and AF segments of random length. At 
an intermediate energy scale (~ J) the spins within each 
segment lock into their local ground state with a variety 
of spin values corresponding to the smallest or largest 
possible total spin for AF or FM segments, respectively. 
Because of the quantum nature of the spins these effec- 
tive degrees of freedom couple rather weakly and behave 
like independent spins for a certain range of tempera- 
ture. They provide the effective spins seen in the Curie 
behavior below the crossover and account for the entropy 
below fceT ~ J/5 0]. At a very low energy scale these 
effective spins will correlate among each other. This pic- 
ture has been confirmed by numerical studies of finite 
size chains which show that the low-energy dynamics is 



described by effective spins of varying sizes coupled by an 
exchange interaction whose strength and sign depend on 
the length and order of the segments ||. This suggests 
that in order to study the low-temperature properties of 
the random spin-^ chain, it is necessary to introduce an 
effective Hamiltonian, 



H — JiSi ■ S 



(1) 



where the strength and sign of Ji as well as the size of 
the spins Si are random. The great generality of the 
Hamiltonian (|l|) suggests that apart from describing the 
low-energy behavior of SraCuPti-^Ir^Og it should be ap- 
plicable to a wider range of disordered ID spin systems. 

We investigate the low-temperature behavior of sys- 
tems described by (Q) by means of a generalized version 
of the real-space renormalization group (RG) method in- 
vented by Dasgupta and Ma (DM) j?J to study random 
bond AF S = | Heisenberg chains (see also |^,||). Let 
us define a link in the chain as two neighboring spins and 
the coupling J connecting them. If a link were isolated 
from the rest of the chain, it would form a local ground 
state of maximum (J < 0) or minimum (J > 0) spin 
with an energy gap A to the first excited multiplct. For 
a FM bond A = -J(S L + Sr) while for an AF bond 
A = J(\Sl — Sr\ + 1), where Sl and Sr are the left 
and right spins of the link, respectively. If we keep track 
of whether a link is ferromagnetic or antiferromagnetic, 
there is a one-to-one correspondence between gaps and 
coupling constants and we characterize from now on a 
link by {A, Sl, Sr}. Now we focus on the strongest link 
in the chain, defined as the link with the largest gap 
A = Ao- If the distribution of gaps is broad, the gaps of 
the two neighboring links, Ai and A2, are typically much 
smaller than Ao so that the two spins Sl and Sr, to a 
good approximation, lock into their local ground state. 
Consequently we replace the link {Ao, Sl, Sr} by a sin- 
gle effective spin S — \Sl ± Sr\ representing the local 
ground state of minimum (AF) or maximum (FM) spin. 
The weaker neighboring bonds are then taken into ac- 
count perturbatively in A12/A0, leading to an effective 
interaction between the spins Si, S, and £2 (see inset in 
Fig. |l|). Spin isotropy is conserved in this procedure, and 
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to first order in Ai^/Ao no next-nearest-neighbor inter- 
actions are generated. The case of S = (corresponding 
to an AF link with Sl = Sr) needs a special treatment, 
since to hrst order this would simply decouple the left 
and right part of the chain. In this case we remove both 
Sl and Sr (which together form a singlet) and make use 
of the fact that second-order processes introduce a weak 
effective coupling between Si and S2. This last step cor- 
responds to the DM treatment of the purely AF spin-i 
model. 

Integrating out the strongest link successively in the 
manner described above preserves the form of the Hamil- 
toninan (Q) but changes the distribution of links and, in 
particular, lowers Ao- Thus the RG procedure generates 
a flow in the distributions of FM and AF links, 

P F < A (A ;A,S L ,^). (2) 

After having integrated out links until the largest remain- 
ing gap in the chain is Ao, the chain is described by the 
Hamiltonian ([!]) with AF and FM links distributed ac- 
cording to (|^) . For low enough energies Ao we expect the 
distributions to flow to fixed-point distributions where 
P F and P A exhibit scaling behavior 

P F ' A (A ; A, S L ,S R ) = ApQ F > A (-*, (3) 

In principle the exponents a, (3, and 7 in P F could dif- 
fer from the ones in P A . However, this is very unlikely, 
because FM and AF links are strongly interrelated with 
each other as we can see from the fact that a large number 
of the spins belong to both an AF and a FM link and that 
in the RG procedure an AF link can be converted into 
a FM link and vice versa. Moreover, different exponents 
7 would imply that one of the probability distributions 
scales to zero, leaving us with a purely FM or AF chain. 
Indeed our numerical RG calculation shows that this is 
not the case. In addition to the exponents in (^|) we are 
interested in how length scales with Ao. Defining n as 
the ratio of the original number of sites to the number of 
effective spins, we expect n ~ A^ s . 

The exponents in (^|) are not all independent. From the 
normalization condition J2s ^s R Jo^-^ + P A )dA = 1 
for any value of Ao follows 7 = /3 — 2a. Furthermore, if 
Q F,A describe fixed-point distributions the average gap 
(A) has to scale linearly in Ao, i.e., (3=1. Finally we 
argue that <5 = 2a. An effective spin at energy scale Ao 
is built up by two more 'primitive' spins which in turn 
are built up by two even more 'primitive' spins etc. The 
total spin of such a unit is just the sum of the participat- 
ing original spins, where each spin enters the sum with 
the same (opposite) sign as its neighbor if the coupling 
is ferrromagnetic (antifcrromagnctic). Thus an effective 
spin made up of n original spins is a sum of n indepen- 
dent random variables and scales as S 2 cx n (see below). 
Thus, we find n«S 2 « Aq 2q and hence S = 2a. 



We have numerically performed the RG for five chains 
with different starting distributions, including one purely 
AF distribution and one distribution simulating the low- 
energy physics of Sr3CuPto.8lro.2O6 Eofl . In each case we 
keep the number of sites fixed to lOMyy adding one link 
after each decimation and iterate up to 2 • 10 7 steps. De- 
tails of this calculation will be presented elsewhere [ pd| . 
In all five chains the distribution after appropriate scal- 
ing eventually approaches the universal fixed-point dis- 
tribution (|J). Integrating Q F and Q A over two variables 
gives the distribution of spins and the distribution of gaps 
plotted in Fig. |l|. The ratio of FM bonds to AF bonds 
is 0.59 at the fixed point and from the scaling of the av- 
erages (see Fig. ||) the exponents are determined to be 
a = 0.22 ± 0.01, (3 = 1.00 ± 0.005 (in both P F and P A ) 
and 5 = 0.44 ± 0.02. Thus the scaling form (|) and the 
relations between the exponents are indeed verified in our 
numerical calculations, and the fixed-point distributions 
become 

pF ' A = in^Q F - A ( A / A °> s ^d, S R A%) (4) 
A o 

and 

n ~ Aq 2q ; a = 0.22 ± 0.01 . (5) 

The scaling forms (^) and ([|) are the key results from 
which we obtain the magnetic susceptibility Xi entropy a, 
and specific heat C. At finite temperature T the RG flows 
will stop at Ao ~ ksT. All pairs of spins in links with the 
gap larger than Ao are frozen to form effective large spins 
according to the RG, while the spins that have survived 
down to Ao are typically much more weakly coupled and 
essentially free. These free spins give a Curie contribu- 
tion to the magnetic susceptibility per unit length, 

X = ^ = cA^ = c 

L 3k B T n TAq 2 " T' { ' 

The Curie constant c can be calculated in the following 
way. Our Hamiltonian (1) contains only nearest-neighbor 
coupling and, therefore, does not give rise to any frus- 
tration. In order to estimate (S 2 ), the average size of 
the effective spins, we consider an isolated system with 
n spins and search for the total spin quantum number 
of its ground state. From the multiplet of the ground 
state we choose the state with the largest z-component, 
S* ot = Stot (total spin size). This state contains the basis 
state corresponding to the classical ground state with all 
spin axes parallel to the z-axis. The correlation of clas- 
sical spins is given by S^S^ = S t Sj ■Y^k= 3 s S n (-^)i 
where Si is the size of the spin i. The disorder averaged 
value of S* ot can now be obtained within the random 
walk picture as, 

n 

((^ ot ) 2 )=^(^ M ) 2 + 2^5 c ^y 

i=0 i>j 
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n Ti—1 i—1 



= n{Sf) + 2{S t ) 2 E E II S 8' n (- J fe), (7) 

i=l j=0 k=j 

where (Si) and (Sf) are averages of the spin size of the 
original spin system (0). Using the average (sgn(— J;)) = 
2q — 1 (q: probability that a bond is ferromagnetic) we 
find the Curie constant, 



c = 



JJ_ 

3k f 



2 r 



2q-l 



{Sif 



(8) 



in the limit T — * (n — ► oo). In particular, in the case 
where all spins are 1/2, the low-temperature Curie con- 
stant is c = (/x 2 /12/cs)g/(l — q), which is in contrast to 
the high-temperature value c = [i 2 /Aks- 

The entropy at zero external magnetic field follows 
from similar arguments. At a given temperature T the 
spins that have not yet frozen into large spins are essen- 
tially free and contribute to the entropy as 

t(T,H = 0) fc B ln(2(S) + l) 2 

= oc (xaksl ml, (9) 

L n 

which give the leading term of the specific heat as 



^p- ocT 2Q |lnT| 
L 



(10) 



In a finite magnetic field H and at finite temperature 
the RG flows are interupted either by the thermal en- 
ergy ksT or by the magnetic energy Eim = (i(S)H. In 
the former case the magnetization is given by (|6|) and 
the entropy by (|^). Otherwise, the magnetic field even- 
tually drives the system away from the fixed point of 
zero magnetic field into a state of aligned effective spins 
where the entropy is zero and the magnetization sat- 
urates. A magnetic field H starts to align the spins 
at an energy scale A ~ fi{S)H, i.e., A ~ H 1 ^ 1+a 
and the saturated magnetization can be estimated to be 
M/L ~ fi{S)/n ~ H a / 1+a . The condition that the chain 
is not yet dominated by thermal fluctuations thus be- 
comes k B T < A ~ and hence 



M (T, H) 
L 



\H/T, 



rpl- 



< bH 
» bH, 



(11) 



where b is a non-universal dimensionful constant. Simi- 
larly the entropy o~(T, H) decreases rapidly towards zero 
at T 1+a ~ fi(Si)H when the magnetic field starts to align 
the spins. 

In the random bond AF spin-i chain the DM RG 
treatment becomes asymptotically exact S because the 
perturbation parameter e = A t /A Q approaches zero in 
the low-energy limit (A 4 : the typical gap in the chain). 
This is not the case in our RG scheme where e ~ 0.2 for 
the fixed-point distribution (see Fig. Q) and higher order 
terms induce interactions other than nearest-neighbor 



coupling. However, such interactions between spins sepa- 
rated by a distance d have couplings of order e d . This ex- 
ponential dependence on d is preserved in the RG so that 
the application of an effective nearest-neighbor Hamil- 
tonian (|l|) is justified. Furthermore, terms of the type 
(Sj • Si+i) m may appear in higher order perturbation. 
Although they may change the local correlations, they 
would not affect the relations between the exponents. 
For this conclusion, it is important to remember that 
the effective spin scales as S oc n 1 / 2 . This scaling results 
from the random correlation of the spins within clusters 
of length n and can be understood immediately within 
a random walk picture, which is quite robust against 
changes of details of the correlation. A consequence of 
this is the relation between the exponents, a = 25 which 
reduces the number of free exponents to one, a. The 
only possible effect of higher order terms is to modify 
slightly the exponent a. In addition, we emphasize that 
the Curie-like behavior of the magnetic susceptibility (||) 
is independent of a. 

The low-energy physics of our system is very differ- 
ent from that of the random AF spin-i Heisenberg chain 
studied in Refs. §-§. In the AF spin- j chain the spins 
successively lock into singlets which afterwards do not 
participate in the low-energy physics. This random sin- 
glet phase |l2| is characterized by the formation of sin- 
glets between spins that may be far apart from each 
other. The coupling between two S = ^ spins that have 
survived down to some low-energy scale is mediated by 
virtually exciting all the singlets in between the two spins 
into triplets. This leads to a coupling that decreases ex- 
ponentially with length, J = A oc exp(— y/n). This, to- 
gether with the fact that the spin size remains constant 
(S = |), leads to susceptibility \ oc T _1 | lnT| -2 and en- 
tropy a(T,H = 0) oc |lnT|~ 2 . In our terminology this 
would correspond to P F = 0, a = 0, and 6 — up to log- 
arithmic corrections. In contrast, our fixed point is char- 
acterized by large spins and a power law relation between 
energy and length leading to (||), (||), (||), and (|Io|). It is 
clear that the two systems belong to different universality 
classes. We emphasize here that the universality class of 
our model includes a large variety of random spin mod- 
els. In particular, an AF spin chain with random spin 
sizes belongs to this class. (It is easy to see that such 
a system generates effective FM interaction among effec- 
tive spins as the RG iteration proceeds.) Similarly, any 
spin chain with randomly distributed AF and FM bonds 
should belong to this universality class |13| . Hence the 
S = 5 AF fixed point studied in Refs. |7]-|9f] is unstable 
toward the presence of FM couplings as well as toward 
randomness in the magnitude of the spins. In either case 
we expect the system to flow to the fixed point discussed 
in this paper. 

An important candidate for this universality class is 
represented by the compound SrsCuPti-^Ir^Oe men- 
tioned earlier, where the low-energy physics is described 
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by the Heisenberg model (|lj). This system shows three 
Curie regimes with different Curie constants At 
high temperature it corresponds to independent S= 1/2 
spins, the intermediate regime is governed by free effec- 
tive spins formed on the FM and AF segments, and the 
low-temperature regime scales to our universality class. 
The slow approach of (n) to the fixed point shown in 
Fig. |2| means that it will be difficult to observe the fixed- 
point value of the specific heat exponent whereas we find 
numerically that the limit of x given in (||) and (||) should 
be more readily observable. 
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FIG. 1. (a) The fixed-point distributions of the spin sizes, 
Pp A {S) = J2 S J dAQ F ' A . In the figure both distributions 

are normalized so that y) g fffi = 1. (b) The fixed-point 
distributions for the gaps, P^' A (A) = ^ s g Q F ' A . The 
distributions are normalized as in (a). Inset: Schematic pic- 
ture of the decimation procedure. 

FIG. 2. Scaling of the averages which determine the ex- 
ponents a, p, and <5. The averages of the gap (A) and spin 
size (S) are calculated from P F (Ao; A, Sl, Sr). The corre- 
sponding averages for P A are almost identical. The initial 
distribution simulates Sr3CuPto.8Ro.2O6 [0. We have found 
more rapid crossover to the fixed point for other initial distri- 
butions, e.g. the fiat distribution P F ' A = const., —1 < J < 1, 
Sl,b e {1/2,1,. .,3}. 
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